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The mechanism of how critical end points of the first-order valence transition (FOVT) are 
controlled by a magnetic field is discussed. We demonstrate that critical temperature is sup- 
pressed to be a quantum critical point (QCP) by a magnetic field. This results explain the 
field dependence of the isostructural FOVT observed in Ce metal and YbInCu4. Magnetic field 
scan can make the system reenter in a critical valence fluctuation region. Even in intermediate- 
valence materials, the QCP is induced by applying a magnetic field, at which magnetic suscep- 
tibility also diverges. The driving force of the field-induced QCP is shown to be a cooperative 
phenomenon of the Zeeman effect and the Kondo effect, which creates a distinct energy scale 
from the Kondo temperature. The key concept is that the closeness to the QCP of the FOVT is 
vital in understanding Ce- and Yb-based heavy- fermions. This explains the peculiar magnetic 
and transport responses in CeYIns (Y=Ir, Rh) and metamagnetic transition in YbXCu4 for 
X=In as well as the sharp contrast between X=Ag and Cd. 
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1. Introduction 

Quantum critical phenomena in itinerant fermion sys- 
tems with strong correlations have attracted much at- 
tention in condensed matter physics. When the contin- 
uous transition temperature of the magnetically ordered 
phase is suppressed by controlling material parameters 
and reaches absolute zero, the quantum critical point 
(QCP) emerges. In the paramagnetic metal phase near 
the QCP, enhanced spin fluctuations cause non- Fermi liq- 
uid behaviors in physical quantities exhibiting quantum 
criticalities, 1-3 and even trigger other instabilities such 
as unconventional superconductivity. So far, the mag- 
netic QCP and the role of spin fluctuations have been 
extensively discussed from both theoretical 1-3 and ex- 
perimental sides. 4 

Recently, critical phenomena associated with charge 
degrees of freedom have attracted attention. In particu- 
lar, valence instability and its critical fluctuations have 
attracted much attention as a possible origin of anoma- 
lies in Ce- and Yb-based heavy-fermion systems. 5,6 Va- 
lence transition phenomena were closely studied four 
decades ago under the label of intermediate valence. Ev- 
idence of its occurrence comes from the 7-a transition of 
Ce metal 7 six decades ago (see Fig. 1(a)) with first-order 
valence transition (FOVT), which in the (temperature 
T, pressure P) plane starts at T v ~ 120 K and termi- 
nates at the critical end point (CEP) (Tcep ~ 600 K, 
P C ep ~ 2 GPa). 8 As the intercept of the T V (P) line of 
the FOVT is rather high at P = 0, no quantum criticality 
was discovered. 9 For many anomalous Ce compounds, no 
FOVT has been detected despite the fact that their va- 
lence deviates from three where the occupation number 



(nf ) of the 4f shell is unity; the conditions are such that 
the system is always in a valence crossover regime, i.e., 
the system escapes from the CEP but, as will be stressed 
later, can feel its proximity. 
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Fig. 1. (color online) Temperature-pressure phase diagram of (a) 
Ce metal 8 and (b) YbInCu4. 13 ' 14 (a) The first-order valence 
transition between the 7-phase and the Q-phase (solid line) ter- 
minates at the critical end point (CEP) (solid circle) in the fee 
lattice. The shaded area around P = represents the /3-phase. 
(b) The first-order valence transition (solid line) is suppressed 
under pressure in the cubic AuBes (C15b type) lattice. Note 
that the T axis is shown on the logarithmic scale. The shaded 
area represents the magnetically ordered phase. 13 rif denotes the 
number of electrons per Ce in (a) and the number of holes per 
Yb in (b) (see text). 



An excellent example of FOVT for Yb systems was 
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reported for YbInCu 4 (see Fig. 1(b)), 10-14 and the Yb 
case can be regarded as the hole analog of Ce; fif being 
the hole occupation number of the 4f shell (14 for Yb 2+ 
and 13 for Yb 3+ where fif = 1). Although YbInCu 4 as 
well as Ce metal is a prototypical example that shows the 
FOVT, most Ce- and Yb-based heavy-fermions seem to 
be in a valence-crossover regime. When the CEP is sup- 
pressed by tuning material parameters and enters the 
Fermi-degeneracy regime, diverging valence fluctuations 
are considered to be coupled with Fermi-surface instabil- 
ity. This multiple instability seems to give the key mech- 
anism that dominates the low-temperature properties of 
materials including valence-fluctuating ions such as Ce 
and Yb. It may play a dominant role in heavy- fermion 
quantum instabilities. 

In Fig. 1(b), with increasing P, T V (P) decreases and 
even becomes sufficiently low such that FOVT reaches 
the magnetic boundary 13 ' 14 clearly in a narrow P win- 
dow. The interplay between valence transition and mag- 
netic transition can be strong. Looking carefully into the 
disappearance under pressure of long range magnetism 
in Ce- and Yb-based heavy-fermion systems, this inter- 
play is often strong. 9 There are only a few material series 
like the CeCu2Si2 15-19 series and CeRu2Si2 9 series where 
the magnetic QCP at P c is not coupled with valence fluc- 
tuation. Thus, discriminating between valence and spin 
quantum criticality is often difficult. 

The proof of valence fluctuations in the quantum de- 
generacy regime seems to be supported by evidence of 
the two-superconductivity mechanism in the CeCu 2 Si2- 
CeCu 2 Ge 2 series where P v - P c ~ 4 GPa. 15-19 A 
marked increase in the superconducting transition tem- 
perature Tsc is observed at a pressure where the valence 
of Ce changes sharply in CeCu 2 Ge2, 15 CeCu2Si2, 16 ' 17 
and CeCu2(Ge x Sii_ x )2- 18 ' 19 The importance of quan- 
tum criticality is shown above Tsc by the observation 
of non- Fermi liquid T-linear resistivity in wide tempera- 
ture region. The T-linear resistivity has been observed in 
a variety of Ce- and Yb-based heavy-fermion systems. 6 

Theoretically, the possibility of the valcnce- 
fluctuation-mediated superconductivity in the P-T 
phase diagram of CeCu2Ge 2 was pointed out in rcf. 5. 
It was shown that the T-linear resistivity emerges 
in a wide temperature range near the QCP of the 
valence transition. 16 Residual resistivity was also shown 
to be markedly enhanced when the system is tuned 
to approach the QCP by controlling P and/or the 
concentration of the chemical doping. 20 Near the QCP, 
the superconducting transition temperature was shown 
to be enhanced by valence fluctuations on the basis 
of the slave-boson mean-field theory taking account of 
its Gaussian fluctuations. 21 The stability and lattice- 
structure dependence of density-fluctuation-mediated 
superconductivity were argued phenomenologically. 22 

Recent numerical studies have revealed the significance 
of valence fluctuations near the QCP of the FOVT 23 and 
clarified its new aspects: The emergence of unconven- 
tional superconductivity due to an anomalous increase 
in the coherence of quasiparticles near the QCP, and the 
absence of phase separation as well as non diverging total 
charge compressibility even at the QCP at least in elec- 



tronic origin due to the non conserving order parameter 
of the valence transition. 23 

In (T, P) phase diagrams of heavy-fermion systems, 
magnetic, valence, and superconductivity boundaries can 
seriously cause interference. This suggests the idea that 
this interplay also occurs in the (T, magnetic field h) 
plane for P close to P v near the QEP. Valence fluctua- 
tions are essentially relative charge fluctuations between 
f and conduction electrons. Hence, it is highly nontrivial 
how magnetic field affects the valence QCP as well as 
QEP. To resolve these fundamental issues, we have the- 
oretically studied the magnetic field dependence of the 
critical points of the FOVT. 24 

In this paper, we report the mechanism of how the 
QCP as well as the CEP of the FOVT is controlled by a 
magnetic field in great detail. We discuss how this newly 
clarified mechanism gives an explanation of unresolved 
observations in Ce- and Yb-based systems. First, we 
show how critical end temperature is suppressed to ab- 
solute zero by applying a magnetic field, which explains 
the field dependence of the isostructural FOVT tempera- 
ture observed in Ce metal and YbInCu 4 . Our results also 
explain the peculiar magnetic response in Celrln 5 , where 
the first-order transition line emerges in the temperature- 
magnetic field phase diagram, giving rise to the increase 
in residual resistivity as well as the appearance of the 
T-linear resistivity. The differences in the location of the 
material with respect to CEP explains the sharp contrast 
between YbAgCu 4 and YbCdCu 4 in their magnetization 
curves in spite of the fact that both have nearly the same 
Kondo temperatures. Our results indicate that the QCP 
as well as the FOVT is induced even in moderately inter- 
mediate valence materials by applying a magnetic field, 
which causes various anomalies such as non-Fermi liq- 
uid behavior in the resistivity, the increase in the resid- 
ual resisitivity, and diverging magnetic susceptibility. We 
discuss the significance of the proximity to the critical 
points of the FOVT to understand unresolved phenom- 
ena in Ce- and Yb-based heavy-fermions. The key con- 
cept is the closeness to the QCP of the FOVT. 

2. First-Order Valence Transition under Mag- 
netic Field at Finite Temperature 

To give a quantitative outlook of the field dependence 
of the valence transition, let us consider the Claudius- 
Clapeyron relation for the FOVT temperature T v : 

dT v _ m K - m M y 
dh Sk — Smv 
where m and 5* denote the magnetization and entropy, 
respectively, and h denotes the magnetic field. Here, K 
indicates the Kondo regime where the f-electron (hole) 
density per site fif is close to 1 in the Ce (Yb)-based 
system, i.e., Ce^^f 1 ) and Yb 3 +(4f 13 ), and MV indi- 
cates the "mixed- valence" regime with fif < 1. Since 
the magnetization as well as the entropy in the Kondo 
regime is larger than those in the MV regime, as ob- 
served in the specific heat and the uniform susceptibil- 
ity, it turns out that T v is suppressed by applying h (see 
Fig. 2(a)). Then, the critical end point is eventually sup- 
pressed to T = K by h. 
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Furthermore, the field dependence of T v in the zero- 
temperature limit is also derived using the above rela- 
tion: For T — ► 0, the entropy shows the T-linear behavior 
in both the Kondo and MV regimes so that 5k — Smv 
is approximated to be proportional to T v in the case 
where T v is smaller than the characteristic energy scales 
in the Kondo and MV regimes. Noting that tok — TiMV is 
temperature-independent for T — > 0, we have ST V /Sh = 



leading to T v = y/2C\\/h v — h with constants 
C\, which explains well the observed behavior in the Ce 
metal 26 and YbInCu4 27 (see Fig. 2(b)). We stress here 
that our analysis not only provides a firm ground for 
small- T v behavior by considering the coherence of elec- 
trons essential for low temperature, but also interpolates 
the high T v satisfying the relation (h/h v ) 2 + (T/T v ) 2 = 
l 26 - 29 to zero temperature, since this relation was derived 
by assuming an isolated atomic entropy, 26 which is jus- 
tified only in the high-temperature regime. Although the 
above discussion is about the FOVT T v temperature, it 
turns out that the critical end temperature Tcep is also 
suppressed, as shown in Figs. 2(a) and 2(b). 



30-36 



(a) T t 




CEP (b) 7i 



QCP 



Fig. 2. (color online) (a) Schematic phase diagram, showing the 
FOVT surface in the T-P-h space, where P represents a con- 
trol parameter (e.g., pressure and chemical concentration). The 
critical end points (CEPs) form a continuous transition line that 
reaches T = as the QCP. (b) FOVT line (h/h v ) 2 + (T/T v ) 2 = 
l26-29 in the T-h plane for a fixed P, corresponding to the dashed 
line in (a). 



According to the Maxwell relation, the volume vari- 
ation V Q (h) with h is related to the pressure derivative 
of the magnetization m (dV a /dh = —dm/dP) propor- 
tional to the Pauli susceptibility in this paramagnetic 
state, and thus inversely proportional to its Kondo tem- 
perature. For Ce, dm/dP decreases under P as fir, thus, 
dV jdh is positive and the system enters in the trivalent 
state upon increasing h. The same occurs for Yb be- 
cause, now, dm/dP increases with P; thus, nf. dV Q /dh 
is negative. The phenomenological dependence of T v (h) 
as reported in Fig. 2(a) was derived in agreement with 
this picture. 29 

3. Extended Periodic Anderson Model 

Although we have shown that the h dependence of the 
critical temperature Xcep can be understood from the 
viewpoint of the free-energy gain induced by the larger 
entropy in the Kondo regime, it is highly nontrivial how 
the QCP of the FOVT is controlled by h at T = 0. To 
proceed with our analysis, we introduce a minimal model 
that describes the essential part of Ce- and Yb-based 



systems 
H = 



2^ EkC^Ck, + Ef 2^ rii 



krx 



N 



N 



Uf^nfa-h^S^ + S?), 



(2) 



i=i 



where Ci a (c| CT ) is the annihilation (creation) operator of 
the conduction electron at the i-th site with a spin <r, and 
fia (ficr) i s that of the f electron. The number operator is 
defined by nf a — a} ia ai a and nf — nf^ + nf^ for a = f and 
c. Here, £k denotes the energy dispersion for conduction 
electrons. Ef is the f level and V is the hybridization 
between the f and conduction electrons. The effect of 
applying pressure is to increase both the hybridization 
V and the f-level Ef relatively to the Fermi level, the 
latter of which plays a more crucial role in the valence 
transition than the former in Ce and Yb compounds. In 
other words, the increase in pressure is parameterized 
essentially by that in Ef. 

The on-site Coulomb repulsion for f electrons is given 
by the U term. The Uf c term is the Coulomb repulsion 
between the f and conduction electrons, which is con- 
sidered to play an important role in the valence tran- 
sition. 21 ' 23 ' 30-36 Namely, the periodic Anderson model 
without Ufc cannot explain a sharp or discontinuous va- 
lence transition as discussed in refs. 21 and 23. For ex- 
ample, in the case of Ce metal that exhibits the y-a tran- 
sition, the 4f- and 5d-electron bands are located at the 
Fermi level. 37 Since both 4f and 5d orbitals are located 
on the same Ce site, this Uf c term cannot be neglected. 
In Yb systems, f ia is regarded as the annihilation 

(creation) operator of the f hole and hence Ef denotes the 
f-hole level. For YbInCu<i, a considerable magnitude of 
the In 5p and Yb 4f hybridization was pointed out by the 
band-structure calculation 38 and recent high-resolution 
photoemission spectra have detected a remarkable in- 
crease in the magnitude of the p-f hybridization at the 
FOVT. 39 These results suggest the importance of both 
Vk and Uf c . 

The reason why the critical-end temperature Tqep is 
so high, that is as much as 600 K, in Ce metal in con- 
trast to that in YbInCu4 can be understood in terms of 
the difference in the magnitude of Uf c . In YbInCu4, Uf c 
originates from the intersite interaction, which should be 
smaller than that of Ce metal. This view also gives the 
reason why most Ce and Yb compounds only show va- 
lence crossover, but not FOVT. Namely, most of the com- 
pounds seem to have a moderate Uf c owing to its intersite 
origin, which is smaller than the critical value for caus- 
ing a discontinuous jump of the valence. However, even 
in the valence-crossover regime, underlying effect of va- 
lence instability causes intriguing phenomena, as shown 
below. It should be noted that the importance of the Uf c 
term in playing a crucial role in the isostructural FOVT 
in YbInCu4 in the hole picture of eq. (2) was pointed 
out in refs. 35 and 36. The last term in eq. (2) is the 
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Zeeman term with h being the magnetic field including 
the g factors. 

3.1 Slave-boson mean field theory 

We apply the slave-boson-mean-field theory 21, 40 to the 
Hamiltonian eq. (2) at T = where the slave-boson- 
mean-field theory is a reasonable approximation. To de- 
scribe the state for U = oo, we consider VfJ a biCi a instead 
of V f\ a a a in eq. (2) by introducing the slave-boson op- 
erator bj at the i-th site to describe the f° state and re- 



quire the constraint 



4bi 



1 by the method 



b]bi - 1 



of the Lagrange multiplier J2i \ n 

For -ffj/tc in eq. (2), we employ mean-field decoupling 
as n\n1 ~ ntnj + n c n\ — ^nfn c . By approximating mean 
fields as uniform ones, i.e., b = (bi) and A = A^, the set 
of mean-field equations is obtained by d(H)/d\ = and 
d(H) jdb = as follows: 



1 

27V 



E 

k(T,± 



1± 



£ftr — £ko 



(4) 



and the following equation holds for the total electron 
number: 



m + n c 



^ £[/(*£) +/(*&)] 



(5) 



Here, f(E) is the Fermi distribution function and E^; 
are the lower (— ) and upper (+) hybridized bands for a 
quasiparticle with spin a, respectively: 



£fa + £k CT ± V (£f cr - £k, 



W 2 



(6) 



where £k<r, £f CT , and V are defined by £k<j = £k + ^fc^f — 
^, E fCT ee ei+\+U lc n c -^f and V" ee 7|6|. The dispersion 
of the conduction electrons is taken as £k = k 2 /(2m) — 
D with —D being set as the bottom of the conduction 
band and the density of states Nq(e) is set to satisfy the 
normalization condition, f^ D deNo(e) = 1 per spin in 
three dimension (see inset of Fig. 3). We take D as the 
energy unit and show the results for V" = 0.5 and the 
total filling n = (fif + n c )/2 = 7/8 below. 

3. 2 Properties at zero magnetic field 

When £f is deep, the Kondo state with fif = 1 is real- 
ized. As £f increases, electrons move from the f level into 
the conduction band via hybridization, giving rise to the 
MV state. Hence, fif decreases gradually for C/ fc = as 
shown in the inset of Fig. 3, as calculated using the slave- 
boson mean-field theory in model (2). As Uf c increases, 
fit decreases sharply as a function of £f. For large U[ c , 
fif = d(H) 1 1 dei shows a discontinuous jump, which in- 
dicates the level crossing of the ground states between 
the Kondo state and the MV state. 21,23 ' 24 (see U {c = 1.6 
and 2.0 in the inset of Fig. 3). The first-order transition 
is caused by Uf c , since a large C/ fc forces the electrons to 
pour into either the f level or the conduction band. 23,24 



2.0 1- a Mixed 

\ Valence 
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Fig. 3. (color online) Ground-state phase diagram in the plane of 
U fc and e f for D = 1, V = 0.5 at n = 7/8. The FOVT line (solid 
line with open triangles) terminates at a QCP (a solid circle) for 
h = 0.00. The dashed line represents the valence-crossover points 
at which Xv has a maximum as a function of Ef for each Uf c . The 
inset shows the energy band of conduction electrons, and fif vs 
e f for U {c = 0.0, 0.4, 0.8, 1.2, 1.6, and 2.0 under h = 0. 



The ground-state phase diagram at a zero magnetic 
field determined by the slave-boson mean-field theory 
is shown in Fig. 3. The FOVT line represented by the 
solid line with open triangles in Fig. 3 satisfies the re- 
lation £f + U{ c n c ~ [i with /i being the chemical po- 
tential 6 ' 21 ' 23 in the mean-field framework. 41 This im- 
plies that the f 1 state with the energy ef + Uf c n c and 
the f° state with a conduction electron at the Fermi 
level with the energy fi are degenerate, giving rise to the 
valence transition. 21 The FOVT line terminates at the 
QCP. The QCP in the £f-C/f c plane is identified to be 
{£? CP , U£ CP ) = (0.356, 1.464), at which the jump in n f 
disappears and the valence susceptibility 



Xv 



d 2 (H) 



dnf 
def 



(7) 



diverges. Namely, valence fluctuations diverge at the 
QCP. Even for Uf c < U® CF , enhanced valence fluctua- 
tions remain, 24 ' 42 as shown by the dashed line in Fig. 3, 
where Xv has a maximum as a function of £f for each Uf c 
(see Fig. 4). The valence-crossover line with enhanced Xv 
regarded as a straight extension of the FOVT line to the 
regime implies that the valence fluctuations 
are a result of the degeneracy of the /° and f 1 states, 
as mentioned above. The characteristic energy scale of 
the system, the so-called Kondo temperature, which is 
defined as 7k = £f<r — Mi is estimated to be 7k = 0.074 
at the QCP. 

Note that the large Fermi surface is realized in both the 
Kondo and MV states: Namely, the number of f electrons 
is always included in the total Fermi volume while keep- 
ing the hybridization between f and conduction electrons, 
as confirmed by the DMRG calculation. 23 This is consis- 
tent with the existence of the QCP in the ground-state 
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Fig. 4. (color online) £f dependence of valence susceptibility Xv 
for U tc from 0.0 to 1.4 at D = 1, V = 0.5, n = 7/8 under h = 0. 
Inset shows Xv vs. £f for U{ c = 1.4. 



phase diagram, since by detouring around the QCP, the 
Kondo and MV states can be adiabatically connected, 
with Luttinger's sum rule satisfied. 

Note that when the hybridization V is decreased (in- 
creased), the QCP is shifted to a larger (smaller) Uf c and 
|ef| position in Fig. 3, as confirmed by the DMRG cal- 
culation as well as the slave-boson mean-field theory. 23 
When pressure is applied to the Ce-based compounds 
such as CeCu2Ge2 and Celrlns, the anion approaches the 
4f electron at the Ce site, which makes the f-electron level 
£f increase. Since the interorbital Coulomb repulsion Uf c 
and the hybridization V also increase under pressure, 
applying pressure is considered to draw a trajectory line 
from the left and bottom position to the right and top 
position in Fig. 3. On the other hand, when pressure is 
applied to Ce metal, it is expected that the increase in 
the hybridization V will be prominent rather than £f, 
because of the monoelemental constitution of the metal. 
Hence, applying pressure is considered to draw a trajec- 
tory line from the left and bottom position to the right 
and top position in the V-Uf c plane, instead of in the 
£f-Uf c plane in Fig. 3. This is consistent with the esti- 
mation of model parameters for the j-a transition in Ce 
metal based on the analysis of photoemission spectra us- 
ing the single-impurity Anderson model. 43 The surface of 
FOVT exists in the parameter space of S{, U{ c , and V for 
the ground state. A trajectory line is drawn in the space 
for the corresponding experimental parameter, such as 
pressure in Fig. 2. 

In Ce metal, the X-ray Lm edge absorption spectra 

led to the conclusion that the Ce valence jumps between 
Ce +3.03 ( 7 phase ) and Ce +3.M ( a phasc ) atT= 3 00 K 44 

One may think that the valence change seems to be too 
small in comparison with our theoretical result shown in 
the inset of Fig. 3. First, we should note that the above 
measurement was performed at a rather high tempera- 
ture (see Fig. 1(a)), and hence the magnitude of the fif 
jump at T = in the inset of Fig. 3 should be markedly 
reduced by thermal fluctuation effects. While, for qual- 
itatively accurate comparison, it is necessary to use re- 



alistic band structures for the f and conduction bands, 
the momentum dependence of the hybridization, and the 
Coulomb interactions Uf c and U in the model (2), it 
should be noted that when the symmetry of the wave 
function of hybridized conduction electrons is the same 
as that of f electrons, the X-ray absorption measurement 
detects the spectra as it comes from f electrons. Hence, 
there is a tendency that this type of measurement un- 
derestimates the magnitude of the valence jump. 

A key parameter for describing valence instabilities 
is the interorbital Coulomb repulsion £/f c , as mentioned 
above. For Ce metal, the onsite Ufc has a considerable 
value giving Ufc > U® CFi5 and hence the FOVT is con- 
sidered to occur at a very high critical end temperature 
Tcep of about 600 K. In the Ce-based compounds such as 
CeCu2Ge2, CeCu2Si2, and CeCu2(Ge x Sii_ x )2, Ufc origi- 
nates from its intersite Coulomb repulsion and hence Ufc 
is reduced from the onsite value, which seems to compa- 
rable to Ufc CP . Namely, these compounds seem to be lo- 
cated in the valence crossover regime (although the sharp 
peak of the residual resistivity and the sharp drop of the 
T 2 coefficient in the resistivity under pressure suggest 
that these are rather close to the QCP 15 ' 16 ' 18 ' 19 ). 

Here, we should also comment on the magnetically or- 
dered phase, which can appear in the ground-state phase 
diagram in Fig. 3 depending on the strength of V. Al- 
though we focus on the nature of the FOVT line with 
the QCP and hence the magnetically ordered phase is 
not shown in Fig. 3, the magnetically ordered phase is 
considered to be realized in the Kondo regime, which is 
basically located in the small-Ef region in Fig. 3. Then, 
in the valence-crossover regime for Ufc < U^ CP , as £{ in- 
creases, the magnetic order is suppressed and the param- 
agnetic metal phase appears. In the paramagnetic metal 
phase, as £f further increases, the Kondo state is changed 
to the MV state at the valence-crossover point repre- 
sented by the dashed line in Fig. 3. We note that in 
the Kondo regime near the QCP in Fig. 3, the super- 
conducting correlation is enhanced, which was shown 
by the slave-boson mean field theory taking into ac- 
count the Gaussian fluctuations 21 and the DMRG cal- 
culation 23 applied to eq. (2). This seems to correspond 
to the T-P phase diagrams of CeCu 2 Ge 2 , 15 CeCu 2 Si 2 , 16 
and CeCu2(Ge x Sii_x)2, 18: 19 where with P application, 
the antiferromagnetic (AF) order is suppressed and in 
the narrow pressure range just before a sharp valence 
increase of Ce, the superconducting transition tempera- 
ture is enhanced. We also note that the reason why the 
superconducting correlation is enhanced was clarified by 
the unbiased calculation: 23 The coherence of electrons 
with large valence fluctuations is enhanced in the Kondo 
regime near the QCP, giving rise to an enhanced pairing 
correlation (see ref. 23 for details). 

We also note that the nontrivial result has been ob- 
tained by the DMRG calculation on the model (2): Total 
charge compressibility is defined by 

^i^V' (8) 

with 2n being the total filling and 2n — fif + n c not 
diverging even at the QCP. 23,42 This is in sharp con- 
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trast to the mean-field result where the phase separation 
is accompanied by the FOVT in the ground-state phase 
diagram, giving rise to diverging n. In the mean field 
framework, the valence fluctuation, i.e., relative charge 
fluctuation diverges at the valence QCP, which triggers 
the total charge instability as well. However, when quan- 
tum fluctuations and electron correlations are taken into 
account correctly, Xv diverges at the valence QCP, but k 
remains finite. 23 ' 42 This has been clarified to be due to 
the fact that the order parameter of the valence tran- 
sition m is not a conserving quantity; [nt,H] ^ 0. 23 
The system can be unstable with respect to the relative 
charge fluctuation while keeping the total charge stable 
(see ref. 23 for details). Hence, it is predicted that, when 
the material parameters could be experimentally tuned 
close to the valence QCP, the compressibility is 

1 (&Vo 



(a) T k 



V \dP 



(9) 



with N e being the total electron number not showing 
divergence at least in electronic origin. 

When we add a temperature axis to Fig. 3 at h = 0, the 
phase diagram of the T-S{-Uf c space is shown schemati- 
cally in Fig. 5(a). The first-order transition surface con- 
tinues to the valence crossover surface. The boundary 
between the two surfaces forms a critical end line, which 
reaches T = K, forming the QCP. At the critical end 
line as well as at the QCP, the valence susceptibility (7) 
diverges, i.e., Xv = °°- Note that, even at the valence- 
crossover surface, valence fluctuations develop well as 
shown in Fig. 4. Figure 5(b) shows a two-dimensional cut 
of Fig. 5(a) for a certain U ic > Ug CP : The FOVT line 
(a solid line) terminates at the critical end point (a filled 
circle), which continues to the valence crossover line (a 
dashed line) . The reason why T v is an increasing function 
of £f can be understood from the Claudius-Clapeyron re- 
lation 



dT v 



(10) 



Sk — Smv 

where n f and S denote the number of f electrons (or f 
holes) per site and the entropy, respectively. Note here 
that the increase in et parameterizes that of pressure in 
our model. We have dT v /def > 0, since n l K > n l MV and 
Sk > Smv are satisfied at least in the deep first-order 
transition region for Uf c > U^ CP . Namely, to achieve 
the free-energy gain caused by the larger entropy, the 
Kondo phase is realized on the higher- T side. Owing to 
the thermodynamic third law, the first-order transition 
temperature T v should be perpendicular to ef for T v — ► 
as understood from eq. (10), which was shown in general 
cases in ref. 46. 

Valence instability is considered to be coupled to vol- 
ume variation. Hence, in our model (2), the effect of the 
hybridization also plays an important role, which might 
share common aspects with the Kondo-volume-collapse 
scenario for Ce metal. 47 However, note that the Kondo- 
volume-collapse scenario assumes a special volume de- 
pendence on the Kondo coupling to realize the first-order 
transition, whose validity should be carefully examined. 




(b) T k 



CEP. "r* 



Fig. 5. (Color online) (a) Schematic phase diagram in the et-U{ c - 
T space for a certain V and U(> Uf c ). The FOVT surface (gray 
surface) with the critical end line (bold line) reaching on T = 
at the QCP continues to the valence-crossover surface (light- 
gray surface), (b) A two-dimensional cut of (a) for a certain 
Uic > ^tc CP - The FOVT line (solid line) T v (e f ) terminates at 
the CEP (filled circle), which continues to the valence-crossover 
line (dashed line) T*(sf). 



Actually, it has been pointed out that the Kondo-volume- 
collapse scenario is not consistent with the isostructural 
FOVT in YbInCu4 by Sarrao: 48 The measured Griineisen 
parameter T = — rflnTK /rfln V = 43 leads to the conclu- 
sion that the difference in Tk at the first-order transition 
should be ATk = 10 K for the measured volume change, 
AVo/Vq = 0.005, which is much less than the measured 
one ATk ~ 400 K (Tk and V denote the Kondo tempera- 
ture and volume, respectively). Thus, the volume change 
is too small to explain the Tk change. In our approach 
(2), the parameters (Uf c ,£{) for each material determine 
whether they show the FOVT, or the valence crossover, 
when T or pressure (or chemical composition) is changed. 
Our scenario does not need to assume the special volume 
dependence on Kondo coupling to cause an FOVT differ- 
ently from that in the Kondo-volume-collapse scenario. 

3. 3 Effects of magnetic field 

3.3.1 Results by slave-boson mean-field theory 

When we apply a magnetic field to the Hamilto- 
nian eq. (2) we find a remarkable result in the valence- 
crossover regime for U lc < U^ CP . Figure 6(a) shows 
the relation of the magnetization m — J2i(Si + 
Sf z )/N vs h for (e f , U ic )= (-0.354, 1.458) (thin line) and 
(—0.349,1.442) (bold line), indicating that the metam- 
agnetism (defined by the diverging magnetic susceptibil- 
ity x = dm/dh — oo) emerges at h = h m = 0.01 and 
0.02, respectively. To clarify its origin, we determine the 
FOVT line as well as the QCP under the magnetic field. 
The result is shown in Fig. 7. It is found that the FOVT 
line extends to the MV regime and the location of the 
QCP shifts to a smaller-J7f c and smaller-|ef| direction, 
when h is applied. This low-h behavior of the FOVT line 
agrees with the low-temperature limit of T v discussed in 
§ 2, in which the FOVT line extends up to the higher 
pressure region as h is increased as shown in Fig. 2(a). 

In Fig. 6(b) we show the m-h curve at U[ c — 1.42 
for £f values ranging from —0.32 to —0.36. The Kondo 
temperature Tk at h = is estimated as 0.0353, 0.0873, 
0.1346, 0.1611, and 0.1823 for e f = -0.36, -0.35, -0.34, 
—0.33, and —0.32, respectively. From these results, the 
mechanism is understood as follows: At h = 0, Tk is 
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Fig. 6. (color online) m-h curve (a) for (e f , Uf c ) = (—0.354, 1.458) 
(thin line) and (—0.349, 1.442) (bold line), and (b) for Ef ranging 
from -0.36 to -0.32 at U ic = 1.42. In both cases, D = 1, V = 0.5 
at n = 7/8. 




-0.37 -0.36 -0.35 -0.34 -0.33 -0.32 



£ f 



originally large for S{ = —0.32 and —0.33, since the sys- 
tem is in the MV regime. However, by applying h, the 
QCP is induced, which makes reduces T K , since the sys- 
tem is forced to be closer to the Kondo regime by h. 
At a magnetic field h — h m where the QCP is reached, 
metamagnetism occurs with the singularity Sm ~ Sh 1 ^ 3 
as shown by Millis et al. 49 (see Fig. 6(a)). The uniform 
susceptibility diverges at the QCP of the FOVT. The 
phenomena lead to the emergence of strong ferromag- 
netic fluctuations. Namely, valence fluctuations diverge 
there, which are essentially charge fluctuations. On the 
other hand, for £f = —0.35 and —0.36, the QCP is not 
reached, so that no metamagnetism appears. 

Note that this mechanism is different from the ordi- 
nary metamagnetism emerging when the magnetic field 
is applied to the Kondo state, which has been dis- 
cussed as the origin of the metamagnetism observed in 
CeRu2Si2. 50-56 Namely, the present metamagnetism is 
caused by the field-induced QCP in the valence-crossover 
regime at h = (for moderate £f and not large Uf c < 
U® CF in Fig. 7) , while the ordinary one is caused in the 
Kondo regime (for deep £f and hence n f ~ 1, see Fig. 7). 
It corresponds to the collapse of antiferromagnetic cor- 
relations and the emergence of ferromagnetic fluctuation 
in a sharp h window. 

An interesting result is shown in Fig. 7, which ex- 
hibits a nonmonotonic h dependence of the QCP: As 
h increases, the QCP shows an upturn at approximately 
h = 0.04, which is comparable to Tk at the QCP for 
h = 0. The upturn of the QCP has also been confirmed 
for a constant density of states No(e) = 1/(21?). This 
nontrivial field dependence of the QCP appears in the 
valence-crossover regime for Uf c < U^ CP in Fig. 7 in 
sharp contrast to the regime for Uf c > U^ CP , where T v 
is monotonically suppressed by h as shown in Fig. 2(b). 

3.3.2 RPA description of QCP 

The nonmonotonic behavior can be understood from 
the structure of the valence susceptibility % v , which is 
given essentially by the RPA, as discussed in ref. 6. 



Fig. 7. (color online) Ground-state phase diagram in the plane 
of Ufc and e f for D = 1 and V = 0.5 at n = 7/8. The FOVT 
line with a QCP for h = 0.00 (open triangle), h = 0.01 (filled 
triangle) h = 0.02 (filled inverse triangle), h = 0.03 (filled star), 
h = 0.04 (filled diamond), h = 0.05 (filled square), and h = 0.06 
(open square). The shaded line connects the QCP under h, which 
is shown as a visual guide. The dashed line represents the valence- 
crossover points at which Xv has a maximum as a function of S{ 
for each Uf c at h = 0.00. 



Namely, it is given as 



xL 0) (*) 



(ii) 



where Xt? ^ s ^ ae bubble diagram composed of f and con- 
duction electrons. In the Kondo regime (h < Tk), where 
f electrons have a predominant spectral weight at approx- 
imately e ~ £f with width A ~ 7rV r2 A(£p), is Mi- 
niated as Xfc"* ~ l/l £ f| an( i i s shown to be an increasing 
function of h. Therefore, U^ CF decreases as h is applied 



QCP, 



QCP 



until it reaches around h ~ Tk, and |e f 
also decreases. For h > Tk, mass enhancement (~ 1/z) 
is quickly suppressed and the MV regime is reached. 

Then, Xt^ ^ s gi ven as Xfc^ ~ 1/A < l/|£f| with us- 
ing the help of shift of the f level towards the Fermi 
level, i.e., £f — ► £f + Uf c Sh c (Sn c being the change in 
the number of conduction electrons per site due to entry 
into the MV regime), so that U® CP becomes larger than 
Ug CP (h~T K ). 

3.3.3 Distinct energy scale from Kondo temperature 

The magnetic field h m at QCP when CEP collapses 
in a magnetic field corresponds to the difference in Tk 
between h = and h = h m : 

h m - Ar£ CP = T$ CP (h ? 0) - T^ CP (h = 0). (12) 

A new energy scale distinct from Tk reproduces the close- 
ness to the valence QCP. Under a magnetic field, the 
proximity of the intermediate- valence crossover regime to 
QCP can lead to the emergence of metamagnetism with 
a jump of m without initially showing the temperature- 
driven FOVT at h = 0. Thus, it is a field-reentrant 
FOVT. 
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Fig. 8. (color online) Magnetization process for ej, = — 2cos(fc), 
V = 1, U = 10 4 , e f = -1 and Uc c = 1 at n = 7/8 calculated by 
the T = DMRG method: (a) m-h curve (filled circle). An arrow 
indicates the metamagnetic transition. Inset: Ef dependence of nf 
extrapolated to the bulk limit for Uf c = (open circle), Uf c = 1 
(filled triangle) and Uf c = 2 (open triangle) at h = 0. An arrow 
indicates Ef = —1. (b) fif (filled triangle) and n c (open square). 
Inset: Lattice structure used in the calculation. In (a) and (b), 
dashed lines represent h = h m . 

3.3.4 DMRG analysis 

To examine the mechanism more precisely, we ap- 
ply the density-matrix-renormalization-group (DMRG) 
method 57 ' 58 to eq. (2) in one dimension. Since valence 
fluctuations are basically ascribed to be of atomic ori- 
gin, the fundamental properties are expected to be cap- 
tured even in one dimension. 23 We show here the re- 
sults for V = 1 and U = 10 4 in eq. (2) at n = 7/8 
on the lattice with N — 40 sites (open-boundary condi- 
tion), as illustrated in the inset of Fig. 8(b). Here, the 
transfer term for conduction electrons is expressed as 
- Etl^L^+i^+H.C.). This lattice may be regarded 
as a one-dimensional mimic of Celrlns and YbXCu4, 
which will be discussed in detail in § 4. 

For h = 0, the relation of fif vs £f for Uf c — 0.0, 1.0, and 
2.0 is shown in the inset of Fig. 8(a). As Uf c increases, the 
change in fif as a function of £f becomes sharp. We show 
in Fig. 8(a) the magnetization m = ^.^Sf 2 + Sf z )/N in 
the MV state for fif indicated by an arrow in the inset of 
Fig. 8(a), which is obtained at h = 0. A plateau appears 
at m = 1 — ri = 1/8, which is expected to disappear if 
we take a more realistic choice of parameters, e.g., the 
momentum dependences of V and £f . The main result is 
that metamagnetism emerges, as indicated by an arrow. 

The increase in fif with a simultaneous decrease in h c 
at h = h m is shown in Fig. 8(b). It is caused by the field- 
induced extension of the QCP to the MV regime. Namely, 
these results indicate that the mean-field conclusion is 
not altered even after properly taking into account the 
quantum fluctuations and electron correlations. 

To further explore the nature of this metamagnetism, 



we calculate T,i( s l z )/ N and E^D/N, and we nnd 
that (Sf z ) decreases slightly at h — h m , while (<Sf z ) 
increases considerably. 24 Since the Kondo cloud is still 
formed even at h = h m , i.e., (S^ • S^) < 0, the decrease of 
(S^ z ) is ascribed to the field-induced Kondo effect, which 
is a consequence of the energy benefit by both the Kondo 
effect and the Zeeman effect. Although this mechanism 
itself has been known to work in the Kondo regime, 55 
this result shows that such a mechanism works in the 
MV regime as a driving force of the field-induced valence 
QCP. 

3. 3. 5 Locality of valence transition 

The DMRG calculation has shown that the magnetic 
susceptibility diverges at the QCP of the FOVT under a 
magnetic field, as shown in Fig. 8(a). It is consistent with 
the slave-boson mean- field theory shown in Fig. 6(a). The 
simultaneous divergence of the magnetic and valence sus- 
ceptibilities at the QCP has been confirmed by the unbi- 
ased calculation; in addition the one-dimensional calcu- 
lation has been shown to be not special. It captures the 
essential physics of valence transition. The main reason is 
the locality of the valence transition. Namely, the valence 
transition has a local atomic origin. Because of this local 
nature, its basic properties and the ground-state phase 
diagram of the valence transition do not depend on spa- 
tial dimensions. Actually, the DMRG calculation in one 
dimension has been known to give essentially the same 
phase diagram determined by the slave boson mean field 
theory (see ref. 23). Recently, the same phase diagram 
has also been obtained by the dynamical mean field the- 
ory in infinite dimension. 59 

The mean-field result supported by the DMRG result 
implies that the RPA approach described in § 3.3.2 is 
qualitatively correct. The validity of the RPA approach is 
also ensured by the perturbation renormalization group 
argument. The dynamical exponent of critical valence 
fluctuations is basically given by Zd = 3, and the con- 
dition d + Zd > 6 is marginally satisfied for d = 3- 
dimcnsional systems. 6 This is a condition for the third- 
order term in the free-energy expansion to be irrelevant, 
as discussed by Hertz for the fourth-order term. 2 Then, 
the universality class of the valence QCP essentially be- 
longs to the Gaussian fixed point, which justifies the RPA 
approach. 

3.3.6 Temperature dependences of FOVT and valence 
crossover under magnetic field 

As shown in § 2, the FOVT temperature T v (h) is 
suppressed by applying h. We note, however, that the 
temperature dependence of T v (h) can change according 
to the location in the ground-state phase diagram (see 
Fig. 7). When the system is located at the deep first-order 
transition side, i.e., at Uf c ^> U^ CP in Fig. 7, T v (h) is a 
decreasing function of h, as shown in Fig. 2(b). Actually, 
the proof by the Claudius-Clapeyron relation is based 
on the fact that ttt-k > timv and <Sk > S'mv at T v (h) in 
eq. (1); this is justified at the deep first-order transition 
side. Ce metal and YbInCu4 correspond to this case. 

On the other hand, near the QCP (Uf c ~ U^ CF ) as 
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well as in the valence-crossover regime at Uf c < U[ 
a different situation can arise: In the case that fif at the 
QCP is not very close to 1 (as is certainly the case for 
Yb systems 60 ), the above relations on magnetization and 
entropy can change. In addition, details of the trajectory 
line of the QCP under h may severely affect the h de- 
pendence of T v (h). For example, Celrln5, which will be 
discussed in detail in the § 4, is considered to be located 
in the valence-crossover regime for U fc < U® CP and the 
field- induced T v (h) is considered to increase under h, in 
contrast to those in the cases of Ce metal and YbInCu4. 

We have already stressed that a valence-crossover sur- 
face exists in the T-ef-U{ c space as shown in Fig. 5(a). 
Hence, even for U( c -C U® CF , the valence crossover sur- 
face will be induced by applying a magnetic field, giving 
rise to an increase in magnetization. The enhancement 
of the magnetic susceptibility at a certain magnetic field 
(see Figs. 7 and 5(a)) will lead to a pseudo-metamagnctic 
effect. Thus, the valence-crossover temperature T*{h) 
also forms a line in the T-h phase diagram. 

When T*(h) is induced by applying h, if the system 
is close to the QCP (namely, is located closely to the 
filled circles in Fig. 7 under h > 0), the following are ex- 
pected to be observed in the T-h phase diagram: (a) The 
T-linear resistivity appears in the wide-T range when h 
approaches h* at which T*(h) becomes zero, (b) Residual 
resistivity is enhanced toward h* and has a maximum at 
h*. (c) Magnetic susceptibility has a peak at T = T*(h). 
(d) NQR frequency changes sharply at T = T*(h), since 
the charge distribution at the Ce or Yb site and its sur- 
rounding ions changes owing to the valence change of Ce 
or Yb, leading to the change in their electric-field gra- 
dient, (e) The lattice constant shows a sharp change at 
T*{h) and hence the magnetostriction changes sharply. 

When the FOVT T v (h) is induced by h, the above 
physical quantities show discontinuous jumps at T = 
T v (h). If the system is close to the QCP, valence- 
fluctuation-induced anomalies such as the T-linear re- 
sistivity will also be observed even in this case. Let us 
conduct a test with experiments. 

In a series of Yb- and Ce-based compounds, the above 
predictions have been actually observed, which will be 
discussed in detail in the next section. 

4. Explanation for YbXCu 4 and CeYIn 5 

We here discuss the potentiality of our theory to re- 
solve outstanding puzzles observed in Yb- and Ce-based 
systems. First, we show how our results explain the 
isostructural FOVT observed in YbInCu4 and the sharp 
contrast between YbAgCu4 and YbCdCu 4 in their mag- 
netic responses. Second, we focus on the peculiar mag- 
netic response in Celrlns, where the first-order transition 
line emerges in the temperature-magnetic-field phase dia- 
gram, giving rise to non-Fermi liquid behavior. Third, we 
argue that the first-order like disappearance of antifcrro- 
magnetism (AF) and the change of de Haas-van Alphcn 
(dHvA) signal observed in CeRhIn 5 at P <~ 2.4 GPa un- 
der a magnetic field h > 10 Tesla may be explained by 
our model. 
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4-1 YbXCu4 systems 

4.1.1 Isostructural FOVT in YbInCu 4 

YbInCu4 is known as a typical Yb compound that ex- 
hibits the isostructural FOVT at T = 42 K 10 - 11 - 61 be- 
tween the high-temperature phase with Yb +2 " 97 and the 
low-temperature phase with Yb +2 ' 84 . 62-64 Namely, in the 
hole picture, fif jumps from 0.97 to 0.84 as temperature 
decreases. This can be understood qualitatively from the 
result shown in Fig. 5(a): in the FOVT region, the larger 
hole-density phase is realized in the high-T phase (the 
Kondo phase), because of the free-energy gain due to 
the larger entropy. Since this high-T phase has a smaller 
f-electron number, the volume of the system is considered 
to be small in comparison with that for the low-T phase. 
Hence, as temperature decreases, the first-order transi- 
tion to the smaller hole-density phase (the MV phase) is 
realized with volume expansion. The collapse of T v under 
a magnetic field was found in macroscopic magnetization 
measurement, explained in phcnomenological approach 
and well confirmed by microscopic X-ray experiments. 64 

As mentioned below eq. (2), the band-structure calcu- 
lations as well as as photoemission measurements suggest 
the importance of the V and Uf c terms in eq. (2) at the 
FOVT in YbInCu4. Here, we should also note the pos- 
sibility that band structures such as semimetallic struc- 
tures also play an important role in the FOVT as pointed 
out in refs. 48 and 65. Although the accurate estimation 
of the values of the model parameters of a model Hamil- 
tonian on the basis of first-principles calculations is an 
important task in the future, we here discuss the basic 
properties of YbXCu 4 on the basis of eq. (2). 

4-1.2 Sharp contrast between YbAgCu^ and YbCdCu^ 

When X=In is used to replace the other elements, 
YbXCu4 docs not show the FOVT, but shows merely 
the valence crossover. For example, X=Ag 66 ~ 68 and 
X=Cd 68 - 69 for YbXCu 4 show neither the FOVT nor the 
magnetic transition and they have the paramagnetic- 
metal ground state. The Kondo temperatures of both 
materials estimated from the magnetic susceptibility 
data 68 are nearly the same, Tk ~ 200 K. A striking 
point is that when the magnetic field is applied to these 
systems, only X=Ag shows a metamagnetic behavior in 
the magnetization curve, while X=Cd merely shows the 
gradual increase in magnetization. 68 

Our results explain why such a sharp increase in mag- 
netization emerges only for X=Ag, but not for X=Cd 
in spite of the fact that both have nearly the same 
Tk's. Figure 9 shows the schematic contour plot of the 
f-hole number per site, fif, which can also be regarded 
as the contour plot of the Kondo temperature Tk in 
the Uf c -6{ plane, because Tk is a function of fif, as in 
Tk oc (1 — fif) /(l — fif /2). 70 In the small-t/f c and small-£f 
region, fif approaches 1, so that Tk becomes small. In the 
large-?7f c and large-£f region, fif is smaller than 1, giving 
rise to a large Tk- As (Uf c , £f) approaches the QCP from 
the valence-crossover regime for U < P , contour lines 
of T K 's get close, and at the FOVT line T K 's show dis- 
continuous jumps. Since the compounds with X=Ag and 
X=Cd have nearly the same Tk, both are considered to 



Full Paper 



10 J. Phys. Soc. Jpn. 




Fig. 9. (Color online) (a) Schematic contour plot of h{, i.e., the 
Kondo temperature Tk in the t/f c -£f plane for the Hamiltonian 
(1) at h = 0. The first-order valence-transition line (solid line) 
terminates at the QCP (filled circle). The arrow represents a 
distance of ~ 40 T between the QCP and YbAgCu 4 . YbCdCu 4 
is located far away from the QCP. (b) Schematic T-E{-U{ c phase 
diagram. YbAgCu4 reaches the valence-crossover surface with a 
distance of about T = 40 K. YbCdCu4 is too far from the QCP 
and hence has too large a distance from the valence-crossover 
surface. 

be located near the same contour area (see Fig. 9). How- 
ever, the compound with X=Ag appears located more 
closely to the QCP of the FOVT than X=Cd, with a 
valence fluctuation energy, i.e., equivalent to a Zeeman 
energy of approximately 40 T. A sharp increase in mag- 
netization will appear in the case of X=Ag; by applying 
h ~ 40 T, the field-induced QCP of the valence transition 
(or sharp valence crossover) is reached. 

A sharp contrast between YbAgCu 4 and YbCdCu 4 
was also observed in the T dependence of the magnetic 
susceptibility x(T)- 68 Although both show nearly the 
same x(0)' s reflecting the fact that both have nearly the 
same Ik's, a broad maximum in x(^) appears at ap- 
proximately T = 40 K only in YbAgCu 4 , but a mono- 
tonic decrease in x(T) appears in YbCdCu 4 as T in- 
creases. This can be understood if YbAgCu 4 reaches 
the valence crossover surface from a distance of approxi- 
mately T = 40 K, while YbCdCu 4 reaches it at too large 
a temperature T interval, as shown in Fig. 9(b). Since the 
magnetic susceptibility x(T) has a peak at the valence 
crossover temperature T* (see Figs. 6(a) and 8(a)), 42 the 
peak of x{T) at T = 40 K in YbAgCu 4 is explained nat- 
urally. Indeed the proximity of YbAgCu 4 to the QCP is 
reflected in the thermal volume expansion directly linked 
to the pressure dependence of the entropy: the volume 
expansion was observed below T = 40 K 71 simultane- 
ously with the sharp valence crossover from 

Y b +2.89 tQ 

Yb+2-87 in YbAgCu 4 , in contrast to YbCdCu 4 . 68 

Hence, the viewpoint of the closeness to the QCP of the 
FOVT expressed in Figs. 9(a) and 9(b) not only explains 
the metamagnetic behavior but also the peak of the uni- 
form susceptibility x{T) consistently. Both phenomena 
are coupled with the local origin of each phenomenon. 

4.2 Celrln 5 

Our results also explain the peculiar magnetic re- 
sponse in Celrln 5 , which shows a jump in the m-h curve 
at 42 T. 72-75 Capan et al. have observed that resid- 
ual resistivity increases, and the Sommerfeld constant in 
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the specific heat shows a diverging increase toward the 
metamagnetic-transition field ~ 25 T. Furthermore, they 
have found that as h increases, the power of the resistiv- 
ity p ~ T a at low temperatures decreases from a = 1.5, 
and seems to approach a = 1.0 judging from the fact that 
the convex curve appears in the p oc T 1,5 plot. 75 Exper- 
imental effort should be made to properly confirm the 
expectation that the T-linear resistivity and the peak of 
the residual resistivity will be observed around h ~ 25 T. 
Capan et al. have pointed out that these anomalous be- 
haviors may be related to the metamagnetic transition 
that forms a first-order-transition line in the T-h phase 
diagram and that this may be the origin of the non Fermi- 
liquid normal state observed at h = 0, although its mech- 
anism has not yet been clarified. 

Our results suggest that the mechanism is the valence 
fluctuation: This can be readily understood if Celrln5 is 
located inside the enclosed area of the QCP line for h ^ 
in Fig. 7. Namely, at h = 0, the system is considered to be 
located in the gradual valence-crossover regime (i.e., for 
U < U® CF in Fig. 7), since no evidence of the first-order 
transition has been observed at any physical quantities 
as a function of T at h = 0. However, when h is applied, 
the QCP of the FOVT reaches and eventually goes across 
the point of the system, causing metamagnetic transition 
in the magnetization curve. Since it has been shown the- 
oretically that the residual resistivity is enhanced near 
the QCP 20 and that the T-linear resistivity is expected 
in a wide-T region, 16 the observed non- Fermi-liquid be- 
havior is quite consistent. Furthermore, the first-order 
transition emerges in the T-h phase diagram in agree- 
ment with our predictions. 

We here remind the readers of the fact that Celrlns 
and CeCoIns have nearly the same crystalline-electric- 
field (CEF) structures 76 and that a change in CEF level 
under a magnetic field cannot explain the metamag- 
netic increase in magnetization in Celrlns, as pointed 
out in ref. 72. Also note that almost the same Fermi 
surfaces in both systems have been obtained by the 
de Haas-van Alphen measurements as well as the first- 
principles band structure calculations. 77 ' 78 However, nei- 
ther the enhancement of residual resistivity nor the 
metamagnetic-transition line in the T-h phase diagram 
has been observed in CeCoIns, in contrast to that in 
the case of Celrlns at zero pressure. These results reem- 
phasize that a distinct energy scale other than the 
Kondo temperature is indispensable for understanding 
the Cell5 systems. 

In order to directly verify the above scenario, detec- 
tion of the Ce valence change at T v (h) in the T-h phase 
diagram is highly desired by measurements such as the 
X-ray adsorption spectra and the NQR electric gradient. 

The notable result is that h scan can lead to a A reen- 
try into the valence critical domain. Our approach as- 
sumes a paramagnetic ground state; in some (J7f c ,£f) 
windows, long-range magnetism will appear as in the 
case of YbInCu 4 under pressure (P > 2 GPa) . As sug- 
gested recently, YbRh 2 Si2 may be a singular spectacular 
case where in the (h, T) phase diagram the magnetism 
and valence fully interact. 60 Such interplay also occurs in 
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CeRhIii5. At zero pressure, CeRhlns is a heavy-fermion 
antiferromagnet with the Neel temperature Tn « 3.5 K. 
At zero magnetic field, at P > 2 GPa, a pure supercon- 
ductivity phase is detected without antiferromagnetism, 
while under a magnetic field reentrant antiferromag- 
netism is detected up to P w 2.4 GPa. 

4.3 CeRhIn 5 

Sharp valence crossover has also been suggested in 
CeRhlns under pressure near P = 2.4 GPa, where the 
resistivity p(T = 2.25 K w T sc ), with T sc being the 
superconducting transition temperature showing a peak, 
as well as the T-linear resistivity in the wide-T range 
emerges. 79-81 Hence, the sharp valence crossover origi- 
nating from the QCP of the valence transition may be 
induced by applying pressure to CeRhlns. Namely, the 
viewpoint of the closeness to the QCP of the FOVT is 
important in elucidating the T-h-P (chemical doping) 
phase diagram of these compounds in a consistent way. 
The sharp peak of p(T = 2.25 K) cannot be explained 
by the many-body correction due to critical AF fluctua- 
tions, 82 but can be understood by the enhanced valence 
fluctuations. 20 

Other lines of evidence for the crucial roles of FOVT 
in CeRhlns under the magnetic field are as follows: 

1) According to rcf. 81, near P = 2.4 GPa under the 
magnetic field of 15 Tesla, both a effective mass of elec- 
trons m* /mo and the coefficient of the T 2 term of the 
resistivity A exhibit a rather sharp enhancement. How- 
ever, m* /mo scales with A 1 / 2 , suggesting that the mass 
enhancement is mainly driven by the "local correlation 
effect" (but not due to critical antiferromagnetic fluctua- 
tions) just as in the case of the metamagnetic transition 
of CeRu2Si2 discussed in refs. 55 and 56. The enhance- 
ment of m*/m can be interpreted as that of the quasi- 
particle density of states near the hybridization gap or 
pseudo gap, which can be approached at approximately 
P = P v ~ 2.4 GPa under the magnetic field as the va- 
lence changes rapidly or discontinuously with pressure. In 
other words, the first-ordcr-like disappearance of antifer- 
romagnetism 81 ' 83 ' 84 and the change of the dHvA signal 
observed in CeRhlns at P = 2.4 GPa under a magnetic 
field larger than 10 Tesla 85 can be naturally understood 
as a FOVT induced by the magnetic field. 

2) According to ref. 81, the upper critical field H c2 ex- 
hibits a rather sharp peak at P = P v where the Fermi 
surface exhibits sharp change from "localized" to "itin- 
erant" under a magnetic field, 85 while the superconduct- 
ing transition temperature T sc is essentially flat around 
P = P v . This fact can be interpreted as an effect of the 
growth of the paring interaction due to the effect of ap- 
proaching the magnetic-ficld-induced critical point of va- 
lence transition. Such behavior reminds us of the case of 
UGe2 at P x = 1.3 GPa in which a magnetic field induces 
the metamagnetic transition between two ferromagnetic 
phases, leading to a sharp increase in H C 2 with a concave 
shape. 86 ' 87 

3) The P dependence of the low-temperature resis- 
tivity p(T = 2.25 K) has a peak at P = 2.4 GPa and 
the emergence of the T-linear dependence of p(T) in the 
vicinity of P = 2.4 GPa 79-81 can be naturally explained 



by the present mechanism. 16 ' 20 

We stress here that the "localized" -to- "itinerant" 
change in electron character reported in the dHvA mea- 
surement 85 can be explained by the Ce-valence jump or 
sharp crossover at P = 2.4 GPa where the number of 
f electrons is always included in the total Fermi vol- 
ume, 23,56 i.e., c-f hybridization is always switched on 
in sharp contrast to that in the Kondo breakdown sce- 
nario. 88 Our mechanism is also consistent with the ex- 
perimental fact that the effective mass of electrons is en- 
hanced even at P — with the Sommerfeld constant 
7 w 56 mJmol -1 K -2 , 79 which is about 10 times en- 
hanced from the LaRhlns value, 78,83 strongly indicating 
the AF state with the c-f hybridization. Furthermore, the 
mass enhancement observed toward P = 2.4 GPa 81,85 in- 
side the AF phase can also be explained by the present 
mechanism. Hence, it should be stressed that the valence 
QCP itself is the source of locality emerging in CcRhIn 5 
without invoking a collapse of Kondo temperature. 

As shown by the phase diagram of CeRh x Iri_ x In 5 , 89 
Cclrln 5 at ambient pressure is moderately far from the 
AF QCP with a distance of about x <~ 0.5. A slight in- 
crease in the nuclear spin-lattice relaxation rate 1/T\T 
indicates that moderate spin fluctuations may exist at 
least at ambient pressure. 89 It has been also reported 
that the magnetotransport measurements under pres- 
sures can be understood from the effects of AF spin fluc- 
tuations. 90 However, the clear difference between Celrln 5 
and CeCoIn 5 emerging in the T-h phase diagram men- 
tioned in § 4.2 cannot be explained only from the sole 
viewpoint of the closeness to the AF QCP. In addition 
to the AF QCP, the influence of the QCP of the FOVT 
is indispensable for the comprehensive understanding. 

Our present viewpoint also gives us a key to resolv- 
ing the outstanding puzzle about the origin of the su- 
perconductivity of Celrln 5 , whose transition tempera- 
ture increases even though AF spin fluctuation is sup- 
pressed under pressure. 89 Since superconductivity will be 
enhanced near the valence QCP, 21,23 the present view- 
point offers a new scenario that the proximity of QCP of 
the FOVT is the main origin of the superconductivity. 
The superconducting window reveals phenomenon other 
than spin fluctuation; the occurrence of superconductiv- 
ity is a unique opportunity for scanning through different 
pairing channels. We have already pointed out in the in- 
troduction that in many heavy fermion compounds even 
for the magnetic QCP the interplay between spin and va- 
lence fluctuations is the main reason for collapse of the 
long-range magnetism. 

4-4 Brief summary 

Detailed discussions for each material have been given 
in from § 4.1 to § 4.3. They are briefly summarized as 
follows: 

1) The field dependences of the FOVT in Ce metal and 
YbInCu4 are clearly explained by our mechanism. The 
effects of the semimatallic band structure on the FOVT 
as well as the qualitative evaluation of Uf c are issues to 
be studied in the future for the complete understanding 
of the valence transition of YbInCu4. 

2) The metamagnetism at h = 40 T and the peak 
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structure in x{ T ) at T = 40 K in YbAgCu 4 but not 
in YbCdCu4 are naturally explained by our mechanism. 
The experimental fact that the valence crossover occurs 
at T = 40 K for h = in YbAgCu4 is consistent with our 
theory. The direct observation of the Yb- valence change 
under a magnetic field at approximately h = 40 T at low 
temperatures is highly desired to confirm our theoretical 
proposal. 

3) We point out that the field-induced FOVT explains 
the T-h phase diagram as well as the non-Fermi-liquid 
critical behavior observed in Cclrln 5 . It has been re- 
ported that magnetotransport measurement under pres- 
sure can be explained by AF spin fluctuations. 90 We 
think that, in addition to the influence of the AF QCP, 
the viewpoint of the closeness to the QCP of the FOVT 
is necessary for the comprehensive understanding of 
Celrln 5 . To examine our theoretical proposal, it is highly 
desired to experimentally determine whether the change 
in Ce valence occurs at the FOVT T v (h) in the h-T phase 
diagram. 

4) We point out that the anomalous behaviors at ap- 
proximately P <~ 2.4 GPa observed in CeRhlns can 
be naturally explained if the FOVT or sharp valence 
crossover of Ce takes place at such a pressure. We think 
that such behavior cannot be explained solely by the AF 
QCP scenario. Our picture gives a natural explanation of 
the origin of the locality as well as the non-Fermi liquid 
behavior without relying on artificial assumptions such 
as the Kondo breakdown. It is highly desired to exper- 
imentally determine whether the Ce valence changes at 
approximately P <~ 2.4 GPa. 

On points 3) and 4), a more quantitative evaluation of 
model parameters including CEF parameters 91 is neces- 
sary for elucidating the P-h-T phase diagram toward a 
complete understanding of the Cell5 system. 

5. Conclusions 

We have clarified the mechanism of novel phenom- 
ena in heavy-fermion systems emerging under a magnetic 
field and have discussed the significance of the proxim- 
ity to the FOVT as a potential origin of the anoma- 
lous electronic properties of Ce- and Yb-based heavy- 
fermions. We have shown that FOVT temperature is 
suppressed by applying a magnetic field, which correctly 
connects the high-temperature result derived from the 
atomic picture of the valence-fluctuating ion to the zero- 
temperature limit consistently with the observations in 
Ce metal and YbInCu4. The important result is that even 
in intermediate-valence materials, by applying a mag- 
netic field, the QCP of the FOVT is induced. The QCP 
shows a nonmonotonic field dependence in the ground- 
state phase diagram, giving rise to the emergence of 
metamagnetism with diverging magnetic susceptibility. 
The driving force of the field-induced QCP is clarified 
to be a cooperative mechanism of the Zeeman effect and 
the Kondo effect, which creates a distinct energy scale 
from the Kondo temperature. 

The use of an extended periodic Anderson model ex- 
plains how quite similar valences may lead to quite differ- 
ent h responses. Our model clarifies why metamagnetic 
behavior appears in YbAgCu 4 but not in YbCdCu4, in 



spite of the fact that both have nearly the same Kondo 
temperatures. The closeness to the QCP of the FOVT 
gives the distinct energy scale, which is a key concept 
to understanding the properties of YbXCu4 (X=In, Ag, 
and Cd) systematically. This viewpoint also explains pe- 
culiar magnetic response in Cclrln 5 where the first-order 
transition line in the T-h phase diagram appears with 
field-induced critical phenomena. The viewpoint of the 
closeness to the QCP of the FOVT is also indispensable 
for understanding CeYIn 5 (Y=In, Co, and Rh) system- 
atically. 

As shown in the present study, the QCP of the FOVT 
and its fluctuations exerts profound influences on Ce- 
and Yb-based materials as a potential origin of anoma- 
lous behavior. Most of such materials are considered to 
be located in the intermediate valence regime, i.e., in the 
region for Uf c < U^ CP in Fig. 7 due to the intersite ori- 
gin of U[ c . However, by applying a magnetic field, the 
valence-crossover surface as well as the critical point is 
induced, which causes various anomalies described in this 
paper. The (P, h) valence transition mechanism clarified 
in this paper can be a key origin of unresolved phenom- 
ena in the family of the materials. 
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